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1. Introduction: 

The Sherrington-Kirkpatrick (SK) model is defined by the Hamiltonian pQ: 

n = --^2j ij a i a j (1) 

where the o~\ are ±1 Ising spins and the couplings Jij are independent Gaussian random 
variables of zero mean and variance equal to 1/N. 

The thermodynamic properties of the model are described by the free energy 
(density) / averaged over the quenched disorder. To overcame the difficulties of 
averaging a logarithm, the average over the disorder is computed using the so called 
replica trick: 

— 7n _ 1 

- PNf = lim (2) 



where (3 = 1/T is the inverse temperature and, as usual, (• • •) denotes the average 
over the disorder. For n integer Z n is the partition functions of n identical, non 
interacting, replicas of the system. The average over disorder couples the different 
replicas. Performing this average, and introducing the auxiliary symmetric replica 
overlap matrix Q ab = jj J^j °~ia,o~ib, with a ^ b, the disorder averaged replicated 
partition functions can be written as pQ: 



^ = Jli\l^ d Q^ NClQ] (3) 

a<b 



with the effective Lagrangian (density): 

= ~^T,^ + m\-n^ (4) 

ab 

fl[Q] = lnTr^ exp ^ Q ab cr Q cr b ^ (5) 

The last term in (jj) follows from the definition Q aa — 1 . The normalization factor in 
([3]) gives a sub-leading contributions for N —¥ oo and is omitted in the following. 

In the thermodynamic limit, N — ¥ oo, the value of the integral in ([3]) is given by 
the stationary point value, and the replica free energy density reads: 

-nf3f = C[Q] (6) 

with Q a b evaluated from the stationary condition 

7 £-£[Q] = 0, a<b (7) 

that is from the self-consistent equation 

Tr CT a a a b exp (^-J2 a b Q ab u a (J b\ 

Qab = -t— v = {<r a (7b), a^b. (8) 

Tr CT exp ( S- Y, ab Qab o a a b \ 

To solve the self-consistent stationary point equation we have to specify the structure 
of the matrix Q a b- This is not straightforward since the symmetry of the replicated 
partition function under replica permutation is broken in the low temperature phase. 
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The Replica Symmetric (RS) Ansatz Q a b = S a b + q (1 — S a b) of Sherrington-Kirkpatrick 
PQ, that assumes the same overlap for any pair of replicas, indeed yields an unphysical 
negative entropy at zero temperature. Following the parameterization introduced 
by Parisi [2] [3J, the overlap matrix Q a b for R breaking in the replica permutation 
symmetry is divided into successive boxes of decreasing size p r , with po = n and 
Pr+1 — 1, along the diagonal, and the elements Q a b of the overlap matrix are assigned 
so that 

Qab = q a nb=r=Qr, T = 0, • • • , R + 1 (9) 

with 1 = Qr+i > Qr > • • • > Qi > Qo- The notation aC\b = r means that a and b 
belong to the same box of size p r but to two distinct boxes of size p r +i < p r - The case 
R = gives back the RS solution, while the opposit limit R — > oo describes a state 
with an infinite, continuum, number of possible spontaneous breaking of the replica 
permutation symmetry. It turns out that a physical solution is obtained only in the 
latter case. Using this structure for Q a b, Parisi and others [2j [3j [4] have shown how to 
obtain solutions with R steps of replica symmetry breaking (RSB) and in particular 
with R — > oo (FRSB), and how to construct equations satisfied by Q{x), the continuous 
limit of the order parameter Q a b for R — > oo 5 . These equations can be solved in 
the full low temperature phase [6j [3 [8j [9] [TUJ [TT]. However working directly with 
Q(x) makes it difficult to keep track, for instance, of the Hessian, and hence of the 
stability of the solution, since the matrix structure of the overlap matrix Q a b is lost in 
the continuous limit The study of the Hessian of the fluctuations around the RSB 
solution with an arbitrary R from the Lagrangean (H])-© is a very hard task. As a 
result stability analysis has been mostly investigated near the critical temperature and 
with the help of a simplified model [TJl [13] , the so called Truncated Model [2j [14] , that 
similarly to the Landau Lagrangian retains only the main mathematical structure of 
the expansion of the replicated free energy in powers of Q a b near T c , where |Q &| <C 1. 

In the present work, we take a different viewpoint and consider the expansion 
of the Lagrangean (HI)-© around the Replica Symmetric ansatz of Sherrington and 
Kirkpatrick. The main motivation for such an expansion is to obtain a simpler 
Lagrangean which, while retaining the replica symmetry breaking properties of the 
original model, is a priori valid in the whole low temperature phase. Anticipating 
our conclusions, we find that the model obtained by truncating the expansion to the 
fourth order, the minimum order required to have a FRSB solution, while improving 
the results obtained from the expansions near T c is valid in a temperature range which 
does not reach zero temperature. 

The outline of the paper is as follows: in Section[2]we construct the approximation 
of il[Q] obtained expanding it around the Replica Symmetric SK solution Q^ S ^ = q 
(a 7^ b) up to fourth order in Q a b — q. The stationarity equation and its solutions 
are discussed in Section [3J The Truncated Model was obtained considering the 
main features of the mathematical structure of the expansion of f2[Q] around the 
paramagnetic solution Q^ 1 * 1 — (a ^ b) to fourth order in Q ab . The parameters 
entering in the model are, however, usually arbitrary and so it is difficult to make 
contact with the original SK model. By using the results of Section[5]we can determine 
the coefficients of the expansion and study the properties of the solution. This is done 
in Section [4] Discussion and conclusions are deferred to Section [5] 

X Paradoxically, it is this continuous limit R — > oo, that imposes the existence of zero modes (at the 
bottom of the replicon bands). Indeed, this limit is necessary to transform the replica permutation 
invariance into a (broken) continuous group thus generating Goldstone zero modes. 
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2. Expansion of the free energy functional around the SK solution: 

To expand the functional Cl[Q] around the SK solution Q a b = q for a =^ b, we consider 
an overlap matrix Q a b of the form 

Qab = S ab + q(l - S ab ) + q ab (10) 
where q is given by the SK Replica Symmetric solution (see below) and q a b the 
deviation from the Replica Symmetric solution. Inserting this form of Q ab into the 
free energy functional (|6]) yields: 

- n PJ = n ~ q - n— q - —q qab - — lab 

ab ab 



lnTr CT exp 



y 5 I 2^ a <* I +yL 9at (Ja(Jb 

\ ab / ab 



+ 0((dl) 



Setting q ab = the above expression leads to the Sherrington-Kirkpatrick free energy 

- /3/sk = y g 2 - + In cosh(/3z) + In 2 + O(n) (12) 

where the overbar denotes the average over the Gaussian variable z: 

+°° dz 



g(z)= I -^=e~* I** g(z). (13) 

Stationarity of /sk with respect to q leads to SK Replica Symmetric solution: 

q = 6 2 ~ : 6> = tanh(/3z). (14) 

For q a b 7^ the free energy functional / can be written, expanding the last term 
in ([TTjl in powers of q a b, as: 

- n/3f = -n/3/sK - y <? H qab ~ If H & 

+Es(i)'((e««)') < («) 

where the subscript "c" indicates that only connected contributions, i.e., only those 
terms that cannot be written as the product of two or more independent sums, must 
be considered. The angular brackets denote the average 



(g(a)) = Y[e^g(cr) + 0(n). (16) 

a=l 

Since <j\ = 1, the last term in (|15p contains only averages of products of spins with 
different replica index. These are easily evaluated yielding 

n h 
a=l 1=1 



[2cosh(j3z)] n ~ h [2sinh(/3z)f + O(n) 



= 6 h + 0{n), ai^---^a h . (17) 

Form the study of the truncated model it is known that terms of order 0(q^ b ) must 
be included into the free energy to break the replica symmetry. Thus in the following 
we shall consider the first four terms of the expansion. 
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2.0.1. Term 0(q ao ): The term of order 0(g a &) is 

( E qab (Ta<Jb / = E 1ab(0- a Crb) = 0' 2 ^ Qab (18) 
\ ab I ab ab 

The choice q — 8 2 , see (Q~4|) , cancels the linear term in the expansion (fPoj) and removes 
the tad-poles. 

2.0.2. Terms 0(q 2 b ): The term of order 0(q 2 b ) reads 

<7ab <7cd {o- a o- b o- c o- d ) (19) 

cd 

To evaluate this term we have to find all different possible ways of equating the ab 
indexes to cd indexes with the constraint, imposed by q aa — 0, that a ^ b and c =/= d. 
There are three possible cases: all indexes different, a pair of equal indexes, and two 
pairs of equal indexes. By noticing that the spin product averages depend only on the 
number of different indexes, and not on the value of the indexes, and that the matrix 
q a b is symmetric, these yield 

2\ 




^2 aa<ib ) ) = 614 E qab qcd + 4 612 E qac qcb + 2 E ? «*» 

\ ab /I abed abc ab 



(20) 



since there are 4 possible ways of equating one index in ab with one index in cd and 
2 was of equating the pair of indexes ab to the pair cd. All sums are restricted to 
different indexes, this is denoted by the prime ""' over the sum sign. Transforming 
the restricted sums into unrestricted ones, i.e., sums over free index, one finally ends 
up with: 

( [^qabO-aO-b] ) = 4 ^ qab qcd + 4 6> 2 (1 - 2 ) ^ gqc gefc 
\ \ ab /I abed abc 

+ 2(1-02)2 (21) 

ab 

This equation has a simple diagrammatic expression. Indeed denoting q a b by a line 
and the vertex where two (or more) indexes are equal by a "dot" , the above equation 
can be written as 

¥ +46» 2 (i -e 2 ) <^ 




+ 2(1-02)2 ^JN, (22 ) 

More details can be found in |Appcndix B| From this form we easily see that the first 
term is a disconnected contribution and hence it does not appears in the free energy 
(fT5|l. therefore to order 0{q a \ b ) the free energy reads 

- npf = -npf SK + q ac qcb + y N E q ab + °( n ^ q3 ab) (23) 

abc ab 

where 

M = 2 2 (1 - 6 2 ), N = (l-0 2 ) 2 -T 2 . (24) 
Notice that the coefficient N is (minus) the Replicon eigenvalue of the Replica 
Symmetric solution [15] . The q a b = solution is hence unstable below T = 1 . 
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2.0.3. Terms 0(q a i b ) and 0(q* b ): These are evaluated as done for the 0(q 2 b ) by 
computing all connected contributions that follows from the expansion of the k = 3 
and k = 4 terms in (115|) . By using a self-explanatory diagrammatic representation 
these are given by: 



where 



P = 246i 2 (l -9 2 ) 2 , Q = 16 6" 4 (1 - 6* 2 ), R = -48 6*2(1 - 6 2 ) 2 , 

J = 16612(1 -6»2)2 ; if = 8 (1 - 6* 2 ) 3 

and 



(26) 
(27) 



\ ab / 



.4 



B 



c <Q- c A +^ d ^- 

3D 0\ + £ V\/ - 2E 0\/ 



with 




+ F 


" + G* {Z^- H 


OO 


(28) 


A = 32 6> 4 (1 


-30 2 )(1- 


e 2 ), 


B = 384 6< 4 (1 - 6*2)2, 


C = 384 6*2(1 - 6*2)3, 


(29) 


D = 64 6*2(1 


— 36» 2 )(1 — 


9 2 ) 2 , 


E = 192 6*2(1 - 6*2)2 


F = 48 (1 - 6* 2 ) 4 , 


(30) 


G = 32(l- 


30 2 ) 2 (l-< 




i? = 96(1- 30 2 )(1- 


6*2)3. 


(31) 



Collecting all contributions up to order 0(g 4 b ), the replica free energy functional reads: 



n f3f = -n(3f SK + 



AT 4 



Mj2<lacqcb + Nj2<l 2 ab 



6(2T 2 )3 



P ^ lacqcdqdb 



'i bed 



abed abc ab abc 

l 



' 24(2T 2 ) 4 

-B ^ q.acqdcq C eq e b + C^] qacqlbQba - C y j q ac qadqd C q. 

abede abc abed 

+4£> ^2 ql c qcb - 3£> ^2 llblbcqbd + E ^2 qabqbcqcdQde 



-A V* 



eQbe QceQde 



cQcdQdb 



abede 



abed 



(32) 



eb 



abc 



abed 



abede 



-2E ^2 qabqbcqld + F ^2 qabqbcqcdqda + G ^ lab - H ^2 lacq 
abed abed ab abc 

-0(n 2 ,q 5 ab ) 
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3. Stationarity equation: 

The equation for q a b follows from the stationarity condition (d/dq a b)f — applied to 
the replica free energy functional (|32|) . In the limit R — > oo this yields 



3(P + Q)S( + R[S 2 + 2Siq(x) + 3Jq{x) 



-6K 



dyq(y)q(y) + s iq (0) 



24(2T 2 ) 4 

b(2Si5 2 - 45? + 2Slq{x)\ + C A(x) 
D(4S 3 - 65i5 2 + l2S iq (x) 2 - 6Sjq(x) 
E(4Sf-4S 1 S 2 -'iSfq(x) 



-4A5| 



-12F 



dyq(y)q(y) 2 + Sfq(0) 



-4Gq(x) 3 -4HS 2 q(x) 



0. 



< x < x c , 



(33) 



where 
A(x) = 2 



dy 3-g(y) 2 q(y) + q 2 (y) ?(») + 5 X g(0) 2 + 5 2 q(0) 



dy 

(4q(x) - 65i 



dyg(j/)<7(2/) + 5ig(0) 



and 



5 n 



dvg(|/)$(tf) 2 -3S? 9 (0)+«(a:) 



da; q(x) n = dx q(x) n - (1 - x c ) q(xj 



(34) 



(35) 



The "dot" indicates the derivative, q(x) — (d/dx)q(x^, while the "hat" the Replica 
Fourier Transform (RFT), that for R —> oo reads [T7] 



q(x) 



q{x) 



dyy^-q{y) - q(x c ), 

dy 
1 ci 

dy- — q(y) + q(0) 
y dy 



RFT 



inverse RFT 



(36) 



(37) 



where q(0) = q(x — 0), and we have neglected the surface term at x = 1 since 
q(x = 1) = q aa = 0. 

Solution of the Stationarity equation 

The complicate integro-diffcrential stationarity equation (|33p can be solved reducing 
it to an ordinary differential equations using the differential operator O = 

§ The RFT was first introduced, directly in the continuum limit (R — > oo) by Mezard and Parisi I18| 
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(l/g(x))(d/dx) to eliminate integrals. Application of O to (|3"3")l leads to 
N 1 



2T 2 3 (2T 2 ) 3 



i?Si + 3 Jq(x) + ZKq{x) 



BS( 



12(2T 2 ) 4 

f[ 1 I dyq(y) q(y) + t?(x) + 4q(x)q{x) - 35 x ?(x) 



2Si g(0) J + D(l2S iq (x) - 35?) - 2ESf + 6Fq{, 



x) 2 



-6Gq(xf - 2HS 2 



= 0. (38) 



The equation is not yet simple enough to be solved. A second application of O, and 
a rearrangement of terms, yields 

8T 2 A(x) + Y(x)q(x) + U(x)q(x) + Z{x)S l = (39) 

where 

X{x) =J-Kx, Y{x) = 2C- AFx, U{x) = AG - 2Cx, Z{x) =AD + Cx (40) 

The integral equation (l39l) can now be transformed into a differential equation dividing 
it by Y(x) and taking the derivative with respect to x. This leads to the first order 
differential equation 

Y(x) [U{x) - Y{x) x] q{x) + (jtq(x) + 8T 2 A + vS x = (41) 

with coefficients 

A = XY — XY = -2CK + AFJ (42) 

/i = UY - UY = -AC 2 + 16FG (43) 

v = ZY - ZY = 2C 2 + 16DF. (44) 

The solution of equation (l4"TT) reads: 

q(x) = T—= J: ~ S a bSi, 0<x<x c , (45) 

y/(X - S) 2 + A 

where 

a = 8T^ b= V - s=^= A=g-, 2 , (46) 
fi n 2b b 

and we have absorbed a factor fi into the definition of the integration constant T. 
The quantity Si is function of T and x c (and temperature), see (|3"5j) . Introducing the 
auxiliary function 

h(z) = f dx X S = + (1 - z) 



y/(x- s y+A y/(z - s) 2 + A 



v/s 2 +A (47) 



_ (z-s)(l-s) + A 
^ - s) 2 + A 
this reads 

Si = Th{x c )-a (48) 

The value of T, and x c , is determined from equations (|38| and (p9| . Replacing 
in equation p9p g(x) with the expression (|45[) yields a linear equation for T. This can 
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be readily solved noticing that since Y does not depend on x we can just set x = 
and use the identity <f(0) = S\. This leads to 

r 



Y\ + Y 2 h{x c ) 



(49) 



where 



a(2G-C-2D), 



Y = AT 2 J{b- 1 

rx=2G(6-l) 7 ^= , (50) 
r 2 = 2Gb - C - 2D. 
Finally the value of x c , for a given the temperature T, is determined from (|38[) . 



Again we can take advantage of the fact that x c does not depend to x and choose 
in a suitable value for x, e.g., x = x c or x — 0. Setting x = into ([55]) a 
straightforward algebra leads to the equation 
1 



2N 



6T 2 



[3 Jq(0) + (R + 3K)Si 



1 



48T 4 

(B-3C-3D- 3E 



(6C+12D)S iq {0) + 

6F)S 2 + (C - 2iT)S 2 ] = 



(51) 



where 

s 2 = - 

and 

h(x) = 



-(1 



(6-1) 



s) 2 + 
-r/i(a;, 



A 



+ / 2 (Zc)-/2(0) + l- 



b{b- 



(6-1) 



^^(x c ) 2 



d.r; 



(x-s) 2 + A 



A tan 



A > 



(52) 



(53) 



Solving equation (|51l) for x c at fixed T yields the value of x c (T), that substituted 
back gives the solution q(x) as function of temperature. In figures [1] and [2] we show 
the solutions Q(x) = q + q{x) for two different temperatures. 

From figures one clearly sees that Q(x = 0) 7^ 0. It grows as the temperature 
decreases, and overcomes q for T < 0.618 . . ., see also figure [3J Retaining in the 
expansion of Q[Q] only terms up to order 0(q^ b ) breaks the replica symmetry, however, 
this approximation is not good enough to change the SK result Q(x — 0) = q 7^ to 
the expected one Q(x — 0) = o(jj| To recover the latter one has to add more terms in 
the expansion, probably all terms. 

Below temperature T = 0.549 . . . equation (|3"5)) ceases to have a physical solution 
and only the SK solution Q{x) = q survives. In figure [3] we show the values of Q(0), 
Q(x c ) and function of temperature. 

I Q(x = 0) must vanish in absence of external fields that break the up/down symmetry. For instance 
in the q > 4 Potts model the symmetry is broken and indeed Q(0) 7^ 
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0.8 



Figure 1. Q(x) versus x at temperature T = 0.9. The horizontal dashed line 
shows the SK solution Q(x) = q. For this temperature we have x c = 0.168846 . . ., 
Q(x B ) = 0.109238 . . . Q(0) = 0.013570 ... and g = 0.102701 . . .. 



Q(x) 




Figure 2. Q(x) versus x at temperature T = 0.7. The full line is the result 
from the expansion around T = 1 to order 0(r 13 ), while the circle are obtained 
from the numerical solution of equation 1 138 1 . The horizontal dashed line shows 
the SK solution Q(x) = q. For this temperature we have x c = 0.3920(6), 
Q(x c ) = 0.3879(1) . . . Q(0) = 0.2232(5) ... and g = 0.3166(5) . . .. 



Sherrington- Kirkpatrick model near T = T c 
1 



11 



0.6 



0.2 



\_ 1 


1 1 1 1 


i 








• • Q(0) 


_ 






■ . Q(x c ) 






w 

\ k 


— q 


- 


- 








- 


\v \ 

X nN\ 






- 
































1 


i , i 


- 

i 





0.2 



0.4 



0.6 



Figure 3. Q(0), Q(X C ) and x c as function of temperature. The full line is the 

SK result q = 8 2 . The Replica Symmetry broken solution ends at temperature 
T = 0.549.... 



3.2. Solution near T c 



1 



Near the critical temperature T c = 1, where both q and q a b vanish, the solution of 
equation (j38]) can be found as a series expansion in the (small) parameter r = T c — T. 
For example to 0(t 5 ) we have 



2r - 4t 2 H 



665 



Q(0) = q 

Q{X C ) = (i 



56 



220 



9(0) = -tt* - — 



68567 
135 7 
3968 



5 + 0(r 6 ) 
- 5 + 0(r 6 ) 



q(x c 



413^ 
"90"' 



O 



r 6 ) 



(54) 
(55) 
(56) 



The resulting series are not convergent, but can be handled by using the Pade 
Approximants. We note that the series expansion of x c has the form of a Stieltjes series 
J2 a n(— T) n . For these series it is known that the diagonal Pade approximant P^(t) 
gives an upper bound and the approximant ^jv+i( t ) a lower bound of the sum |20l . 
Moreover in the limit of large N both approximants converge, and if they converge to 
the same limit this is the value of the sum. By using the Pade approximants we were 
able to use the series expansion almost everywhere in the low temperature phase, 
where the Replica Symmetry broken solution exists. For example for temperature 
T = 0.7 by using the series expansion to 0(r 13 ) we have x c — 0.3920(6), the error 
being estimated from the difference between the Pade Approximants P$ and P$ +1 - 
A comparison between the numerical and the power series solutions is shown in figure 
[U the agreement is rather good. 
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4. Expansion around the Paramagnetic Solution: 



At the critical point T c = 1 the order parameter function Q(x) vanishes, one can 
then think of expanding the functional il[Q] around Q a b = 0, i.e., the paramagnetic 
solution. Such an expansion, first considered by Bray and Moore []3] , is at the basis of 
the so-called Truncated Model [2] largely used to study the properties of the solution 
Q{x) near the critical point. See [19] for an extension to more general models. Despite 
it usefulness, the Truncated Model is a poor approximation for the SK model. Indeed, 
in the same spirit of the Landau Theory of second order transition, it retains only the 
main mathematical structure of the order 0(Q\ b ) expansion of ft[Q] around Qafc = 0, 
but with arbitrary coefficients. Using the results of Section^ we can investigate the 
properties of the 0(Q\ b ) approximation of the SK model. 

The expansion of the replica free energy functional around Q a b — is obtained 
by setting q = 9 = in flU). This yields 

710 f = 71 In 2 + ^ N ^ q 2 ab + + K QacQcbOba 



24(2T 2 ) 4 



6(2T 2 ) 3 

F ^2 labqbcqcdqda + G qtb ~ H q 2 ac q 2 c 



abed 



ab 



(57) 



with 



N = l-T 2 , K = 8, F = 48, G = 32, H = 96 (58) 

Stationarity of (|57[) with respect to variations of q a b leads to the stationary point 
equation, that for R — > oo reads: 

2Nq ^ + ^(J dyq(y)q(y) + S iq (0^j 

dyq(y)q(y) 2 + S 2 q(0) 



24T 4 



3F 



-Gq(x) 3 - HS 2 q(x) 



0. 



(59) 



Applying the differential operator O = (l/q(x))(d/dx), as done in Section [3TT1 reduces 
the above equation to: 



3Fq(x) 2 + 3Gq(x) 2 - HS 2 



0. 



(60) 



This equation is not yet simple enough to be solved. A second application of O leads 
to 



2T 2 Kx - Fxq(x) + Gq(x) = 0. 



(61) 



Dividing this equation by Fx and taking the derivative with respect to x transform 
the integral equation (|61|) into a differential equation, that solved yields: 

q(x) = r = — , < x < x c (62) 



with 



r = it 



2 K y/x* + G/F 
F x c + G/F 



(63) 
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determined inserting the form (lB"2"j) of q(x) into (|6"Tj) . 

The end point x c is not a free parameter and must be determined as function of 
temperature from (|60[) . Inserting from (|62[) and 



- 2T^ 



a; F 

into (|60|) one ends up with the following equation 
K 2 K 2 x 2 + G/F 



(64) 



2N- 



2F 6F (x c + G/F) 2 
Hi,- ' G 



3G 



tan 




= 



(65) 



that solved for x c gives the value of x c (T). 

At the critical temperature T = T c = 1, where N = 0, x c vanishes and increases 
as the temperature is decreased below T c . Introducing the small parameter r = 1 — T, 
the solution of equation (|65[) can be expressed as a power series. For example to order 
0(t 5 ) we have, 



and 



x c = 2r + 12t 2 
Q{x c ) =t + t 2 



v 1 5 

r = —= — -=t- 

a/6 V6 



280 



44 



2437 



701 



37641 



+ 0(r 6 ) 



30763 
30 



r 5 + 0(r 6 ) 



(66) 
(67) 



_L T 2 _ ~3 

Vl3 VG 



2972 



-0(r 6 



Notice that in this case Q(x = 0) = 0. 

The maximum allowed value of x c is 1. Setting x c = 1 into (|65|) . and replacing 
the constants by their values , we find that the Replica Symmetry Broken solution 
(|62p becomes non-physical below the temperature 



FRSB 



21 - V6 tan" 1 




0.9148. 



where x c > 1. At this temperature Q(x c ) reaches its maximum value 



lim Q(x c ) = 

T— >t FR 3 B 225 



21 - \/6 tan -1 




0.16737. 



(69) 



(70) 



We conclude this Section noticing that for temperatures above, but close to, Tfrsb 
equation (|65p can be solved as power series of 1 — x c . We do not report the expansion 
here. 



5. Discussion and Conclusions: 



In this work, we have derived the expansion of the Shcrrington-Kirkpatrick model 
replica free energy functional around the Replica Symmetric (RS) solution Q^ S ^ = 
S a b + <z(l — 5 ao ). We have considered in detail the approximation obtained by 
truncating the expansion to fourth order in Q a i, — Q ab , i.e., the lowest nontrivial 
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approximation to have a continuous Replica Symmetry Breaking. The stationarity 
equation (j3"3")l associated with the approximate free energy functional (j3"2l can be 
solved and the explicit form of the Full Replica Symmetry Broken (FRSB) solution 
Q(x), for < x < be determined. The FRSB solution appears at the critical 

temperature T c = 1, as the RS solution, and exists only down to the finite temperature 
T = 0.549 . . .. Below only the RS solution survives. 

A peculiar feature of the FRSB solution is that Q(x = 0) 7^ 0, and vanishes as 
(T c — T) 3 as the temperature T approaches the critical temperature T c = 1. This 
property can be traced back to the fact that the FRSB solution "opens" around the 
RS solution Q( RS >(x), i.e., Q(x = 0) < q < Q(x c ), as the temperature decreases below 
the critical temperature T c . As the temperature is decreased below T c the RS solution 
q increases and drags Q{x = 0) to finite values. We note that at T = 0.618 .. . the 
value of Q{x = 0) eventually overcomes that of q. 

Setting q — one recovers the expansion of the replica free energy functional 
around the paramagnetic solution = 5 a b to order 0{Q a \ b ). This turns out 

to be rather interesting because such an expansion is at the basis of the Truncated 
Model used to study the properties of the FRSB solution Q(x) near the transition. 
To our knowledge, a study of this approximation with the correct coefficients of the 
expansion was never done. Indeed the Truncated Model, and the one in which one 
keeps all terms generated by the expansion of Q,[Q] to order 0(Q ab ), have been studied 
with arbitrary coefficient. As consequence of this the existence od a FRSB solution 
was always taken for granted, but never verified. We have studied the existence of the 
FRSB solution for this expansion in the last part of this work. Surprisingly it turns 
out that the FRSB solution exists only close to the critical temperature, in the range 
of temperature 0.9148 . . . < T < 1. Therefore such expansions truncated to the forth 
order cannot be used to study the solution of the SK model near zero temperature. 

To summarize, 

• The expansion (to 4 th order) around the (Replica Symmetric) SK solution 

Q^ab^ ~ ^ ab + 3(1 — $ ab ) yields, in the limit R — >• 00, a Q(x) that does not 
vanish for x null, in contrast with the exact Parisi Solution. The solution exists 
only in the range 0.549 . . . < T < T c = 1. 

• The expansion (to 4 th order) around the paramagnetic solution q — 0, yields 
a Q(x) that does vanish for x null. But it exists only close to T c = 1, 
0.915... < T< T c = 1. 

In this work we have studied the existence of FRSB solutions. Finite RSB 
solutions may also exist. These, however, may exhibit problems similar to those 
found for the FRSB solution. For example, inserting the Replica Symmetric Ansatz 
Qab = Q (1 — S a b) into the free energy functional (|57j) . and taking the limit n — ¥ 0, or 
expanding the SK free energy (|12p around q = to the fourth order in q, one ends up 
with 

- PS ^ln2- l -^q^ + ± E q^-^ fE q\ (71) 

Stationarity respect to variations of q yields the paramagnetic solution q — and the 
RS solution 

T 2 - 



q= 17 



3 - ^3(17T 2 - 14) . (72) 
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The latter correctly vanishes at the critical point T = T c = 1, but exists only down 
to temperature T = y/U/17 ~ 0.907 . . ., slightly below the lower end T = 0.915 . . . of 
the FRSB solutions, where the quantity under the square root becomes negative. 

The situation is only slightly better considering the expansion around the Replica 
Symmetric SK solution since now the RS solution exists down to T = 0. The RS 
Ansatz q ao = Sq (1 — S ao ) yields indeed, besides the trivial solution Sq — 0, a Sq ^ 
solution leading to the unphysical result Q = q + Sq ~ — ^(1 — T) 3 as T ~ 1 _ . The 
reason is that the expansion around the SK solution includes the contribution of more 
diagrams: all diagrams that are needed to build the SK free energy. In this sense this 
expansion is a better approximation for the SK model, as also reflected by the larger 
temperature range were the FRSB solution exist. 

The improvement is only apparent since for both expansions the RS solution, as 
well as the Paramagnetic solution for T < 1, has a negative Replicon mass and is, 
hence, unstable. For what concerns the FRSB solution, whether in the full expansion 
(|57p or in the truncated model, they both have the same stability properties in the 
most dangerous sector, i.e., in the Replicon subspace (by virtue of the Ward-Takahashi 
identities [TB]). Thus the FRSB, where it does exist, is marginally stable with null 
Replicon masses. We believe that this feature remains true for the expansion around 
the SK solution as well. 

Despite this, the limited range of temperature (and not including T — 0) in which 
these expansions to 4 th order exist, makes them of little help to study the properties of 
the SK model near zero temperature. To extend the range of validity one should retain 
more terms in the expansion, and probably all terms (or infinite subseries thereof) 
since the particular structures of the expansion (in powers of 0) may otherwise lead 
to difficulties for very low temperatures. We observe that to overcome this problem a 
construction based upon an expansion around a spherical approximation, which leads 
instead to an expansion in T, has been recently proposed |21j . 
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Appendix A. The Truncated Model 



The Truncated model is defined by the free energy 

ab abc ab 

with t = 1 — T and w and u arbitrary and positive. Comparison of (|A.1[) and (|57p 
shows that 

N K - G 

T - 2 T3' W "8T5' U ~32T7 (A - 2) 

while F = H = 0. We can then read the equation for q(x) directly from Section 
[U Setting F = in (|6T|) one readily obtains the known linear form of q(x) for the 
truncated model: 

q{x) = 2T 2 ^x = ^-x, 0<x c <x c (A.3) 
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Finally setting F = H = in (foX))) . and using the above linear form of q[x), yields 
K 2 K 2 

2iV-— x c + — x^O (A.4) 

that gives x c as function of temperature. The value of x c is zero for T = 1 and 
increases as T decreases below 1. By setting x c = 1 into (|A.4jl leads to the critical 
temperature 



if 2 

^- = V 1 -4G (A - 5) 
below which the Replica Symmetry broken solution ceases to exist. If 1 — K 2 /AG < 
the solution exists down to T = 0. If we use the values K = 8 and G = 32 we have 
T trm = 1/V2 = 0.707.... 

We note that due to the presence of T-factors in the relation between (N, K, G) 
and (r, w, u) the critical temperature has a slightly different form if expressed in the 
latter: 

w 2 

T trm = 1 — -7— (A.6) 

4u 

and is valid if w and u are temperature independent. 



Appendix B. Terms 0(q 2 b ): details 
The terms of order 0(q 2 b ) are given by 




(B.l) 



cd 



To evaluate this term we have to find all possible ways of equating the ab indexes to cd 
indexes, with the constraint a ^ b and c ^ d since q aa = 0. In the following to denote 
that a group of indexes must be all different we shall write them in parenthesis, hence 
in the present case we have two write (ab) and (cd). 

We clearly have three possible cases: all different, one equal, two equals. For later 
use it is useful to represent them graphically. If we denote q a t, by a straight line then 
the case of all different indexes is represented as 



c dl _ (B-2) 
and the value of the average is 8 4 since all spin indexes are different. Next there are 
four possible ways of equating one (ab) index to one (cd) index. These are 



c d I \ c dl \ c dl \ c dl (^-3) 

where the indexes connected by a dashed line are equal. In this case two indexes in 
the average are equal, so only two spins survive and the spin average gives 9 2 . 
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Finally there are two possible way of equating indexes (ab) and indexes (cd) with 
the constraint a ^ b and c 7^ d and reads 



c d I \ c d I (B-4) 

In this cases we have two pairs of equal indexes in the average, so all spins disappears 
and the average give 1. 

To evaluate (|B.1[) we have to sum each diagram over a, b, c, d, then it is easy to 
realize that since the matrix q ao is symmetric all four diagrams in (|B.3|) give the same 
contribution, and so do the two diagrams in (IB.4[) . These will be denoted as 




(B.5) 

respectively. 

Collecting all terms we have 



^ q a b o-aVb ) = A g a b q c d + 4 ( 

\ ab J I (abed) 



2 5> QC <z c6 + 2^ b (B 6) 

(abed) (abc) ( a b) 

The restricted sums can be transformed into unrestricted sums by inserting a factor 
(1 — S a b) for each pair of indexes (ab) to enforce the constraint and removing the 
constraint over the indexes. By expanding now the resulting products of (1 — <5)'s 
each restricted sum is finally expressed as a combination of unrestricted sums. 
Diagrammatically we have 




- _ ~ 4 < +2 o 

(B.7) 
and 




(B. 

Inserting these expressions into (|B.6j) after simple manipulations we end up with 

~~ " " (B.9) 



(^Qab(r a a^j ^=9* + A9^(l-9 2 ) <^+2(l-9^ 



The first term is disconnected and hence it does not contribute to the free energy, 
therefore to order 0(q 2 b ) the free energy reads 

- nPf = -n/3/sK - ^ [l ~ /? 2 (1-^ 2 ) 2 ] ]T q\ b 

ab 

+y« 2 (l-« 2 )Ewc6 + 0(^) (B.10) 
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